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^ ■ Abstract 

A theory of evanescent mode atomic mirrors utilising a metallic sheet on a dielectric 

>' 
O ' substrate is described. The emphasis here is on the role of the metallic sheet and on 

m 

O 

o 

O ' intensity, the atomic reflection process is controlled by two separate mechanisms both 

O' 

c/3 \ of which are modified by the presence of the metallic sheet and influenced by the use of 

• 1—1 
^! the field-dipole orientation picture. The first mechanism involves the spontaneous force 

^ : 

p |. which accelerates the atom parallel to the sheet plane. The second mechanism involves 

J . the combined dipole plus van der Waals force which acts to repel the atom from the 

>< 
j^ ■ surface, decelerating its motion until it attains an instantaneous halt before changing 

direction away from the surface at an appropriate turning point in the trajectory. 

Various quantitative features arising from varying the controllable parameters of the 

system, including screening effects as well as desirable enhancement effects, are pointed 

out and discussed. 
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One of the primary aims, yet to be fully accomplished in the rapidly developing field 
of atom optics is the standardisation of the atomic mirror as a routine atom optical 
element. The idea of an atomic mirror, a device capable of reflecting atoms, was 
first put forward by Cooke and Hill in 1982 [1]. Since then, a number of studies, 
both theoretical and experimental, have sought to explore various aspects of atomic 
mirrors, based on different mechanisms of atom reflection from planar surfaces [2-19] . In 
particular, the evanescent mode atomic mirror is suitable for neutral atoms undergoing 
electric dipole transitions [14,18,19]. The mechanism in this type of mirror makes use 
of light evanescing into the vacuum region outside the planar surface of a dielectric 
when laser light is internally reflected. The evanescent light sets up a repulsive dipole 
potential which acts on any neutral atom possessing a transition frequency at near 
resonance with, and blue-detuned from, the laser frequency. The limiting factors of 
such a mirror stem primarily from heating effects and also from fluctuations due to the 
light which become negligible at low intensities. A device operating at low intensities 
was thus needed. This was the motivation behind the studies which highlight the role 
of a thin metallic layer deposited on the surface in providing an enhancement of the 
evanescent component [18,19]. 

The purpose of this article is to report a theory of evanescent mode atomic mirrors 
with a metallic sheet, exploring the nature of the enhancement, determining the general 
features of the mirror and evaluating the trajectories for a typical set of parameters. 

It turns out that a quantitative theoretical analysis demands the implementation of a 
number of steps. First, a mode normalisation procedure in terms of intensity is needed 
for the incident light mode responsible for the generation of the evanescent component. 
Secondly, the jump conditions should be implemented involving the surface current 
arising from the flnite two-dimensional conductivity of the metallic sheet. Thirdly, the 
attractive force between the atom and the dielectric in the presence of the metallic 
sheet should be included in the dynamics to determine the trajectories. Finally, use 



should made of the field-dipole orientation picture for evaluating the radiation forces 
as well as the van der Waals-type atom-surface force. The role of the field-dipole 
orientation picture is to determine an average value for the local orientation of the 
electric dipole moment vector in the presence of the evanescent field. We show here 
that a programme incorporating the above theoretical features permits the atomic 
trajectories to be determined by direct solution of the equation of motion, leading either 
to a reflection of the atom off the mirror or a collision with it, in a manner dependent 
on the chosen set of parameters. Existing treatments determining trajectories are 
primarily based on Monte Carlo techniques. 

The basic elements comprising the atomic mirror are shown schematically in Fig. 1. 
Here a metallic sheet in the form of an inflnitesimally thin layer is deposited on the 
planar surface of a dielectric substrate (or glass prism). Light of frequency uj incident 
from within the dielectric is internally reflected at the interface between the dielectric 
substrate and the metallic sheet. This creates a fleld in the vacuum region which is 
decaying with distance from the metallic sheet and is propagating along the surface. 
A neutral atom possessing a transition frequency uoq < uo and which is moving in 
the plane of incidence would be subject to the repulsive dipole force plus an attractive 
atom-surface force and will also experience a light pressure force parallel to the surface. 
As we show below the combined influence of these forces can be made to control the 
reflection process. 

The relevant electric fleld vector at frequency uo can be written in terms of incident (I), 
reflected (R) and evanescent (1) parts as follows 

E(, r, t) = {(E,(k||, r, t) + E^(ki,, r, t)) e{-z) + Ei(k||, r, t)e{z)} a + h.c. (1) 

where 6 is the unit step function and the flelds are given by 

Ej(k||,r,t) = Aj ( 1,0,-- ^ ) e^^-^vC^ii.rii-o.i) (3) 

V k^2 J 



ER(k|,, r, t) =Ar(i, 0, -^"j e-'^-^e^C^Kii-*) (3) 

Ei(k||,r,t) = B ( 1,0,^ 1 e-'=-iV(''ii-'"ii-^*) (4) 

Here ky is the wavevector parallel to the surface. Its magnitude k\\ is given by c^fcy = 
uj'^e2 sin^ (p where is the angle of incidence. The three quantities between the brackets 
in each of Eqs. @) to @) stand for the vector components parallel to ky , perpendicular 
to it on the surface plane and along the z-direction, respectively, kzi and kz2 (both 
real) are defined by 

kl, = kl-'-^>Q- kl,= '-^-kl>0 (5) 

II ^Z gZ II 

Finally, Aj, An and B are field amplitude factors, to be determined. The notation 
is such that parameters associated with the substrate are labelled by the subscript 
2, while for the outer region (vacuum) the label is 1. Both dielectric functions ei 
and €2 are assumed to be frequency-independent and we take ei = 1, as appropriate 
for vacuum. The position vector is written as r = (r||,2;) in terms of an in-plane 
position vector ry and a z-coordinate relative to the metallic sheet. The role of the 
metallic sheet is primarily to provide a two-dimensional charge density Us and, so, an 
electrical conductivity ingC^ /m*{uj + 27) where m* and e are the electronic effective 
mass and charge and 7^0; accounts for metallic plasma loss effects. The metallic 
sheet only enters the formalism via the electromagnetic boundary conditions involving 
the tangential components of the magnetic fields corresponding to Eqs.(|^) to (^, which 
can be calculated using Maxwell's equation H = —{ieoC^/uj)'Vx E. Application of the 
first boundary condition, namely the continuity of the tangential component of the 
electric field vector at z = 0, yields 

Ai + Ar = B (6) 

The second electromagnetic boundary condition is that the tangential component of 
the magnetic field vector experiences a discontinuity at z = arising from the surface 



current induced by the in-plane component of the electric field at the metallic sheet. 
We have 

2 

if||(0_) - /7||(0+) = , ' . . EM) (7) 

where 0± are the limits as .^ — ^ of (0 ± ^) . Application of this boundary condition 
leads to a second relation connecting the field amplitudes 

-— {Aj - Ar) + t-;—B 



kz2 kzi eom*iu{iu + i^y) 

Elimination of Aji between Eqs.® and (|^) yields straightforwardly 

kz2 ( n^e^ ei 



B = 2Ai 



+ 1 



(9) 



62 \eom*uj{uj + i'y) k^i 
The next step is to determine the value for the amplitude Aj. This is the amplitude 
of the incident field in the unbounded bulk of material 2 and the value of Ai is such 
that the Hamiltonian Tii reduces to the canonical form 



Hi = — dr < e2E] + - — —Hj ^ = -huo [aa) + a) a 



(10) 



Thus we find 
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where V is the (large) volume of material 2. 



The average radiation force acting on an atom of transition frequency ujq moving in 
the vacuum region at velocity v is given by the well known expression [20,22] 

F r, V = 2/1 <^ — -^ 1 -^^=F, + Fd 12 

where F^ corresponds to the first term, identified as the spontaneous force along the 
wave propagation direction ky, and F^ corresponds to the second term, identified as 
the dipole force; A is the dynamic detuning given by 

A = Ao-kj|.v (13) 

with Aq = u — ujq the static detuning of the light from the atomic resonance. The 
force given in Eq.(|l2D is applicable in the low intensity limit for which the saturation 



parameter S satisfies the inequality S = 2f2|./(A^ + F^) < 1, where r2/j(0) is the Rabi 
frequency, defined below, evaluated at 2; = 0, and T is the spontaneous emission rate. 
At the high intensity regime, corresponding to 5 >> 1 one needs to adopt the dressed 
atom approach [21]. The low intensity regime is the correct regime in the context of 
the atomic mirrors considered here for which it can readily be verified that S* < 1, so 
that the dynamics based on the average radiation forces given by Eq. (|T2p is applicable. 

Of the two average radiation forces, it is easy to see that F^ can act as a repulsive force 
provided that the detuning Aq is positive (blue detuning). Note that v, the velocity 
of the atom, is in the plane of incidence as in Fig. 1, i.e. has only two components, a 
z-component and a component parallel to ky. The Rabi frequency Q^ entering Eq.(|T^ 
is defined for the electric dipole /j, in the evanescent field Ei(k||,r, t) 

a/^.Ei(k||, r, t) 



n 



R 



(14) 
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where a is a complex amplitude factor such that a ^ a in the classical electromagnetic 
field limit. It is in fact related to the intensity / of the incident beam by the well known 
relation [23] 

It can be seen from Eq.(4) that the evanescent field possesses two vector components. 
However, once the evanescent field has been set up, the average atomic dipole moment 
vector at any given point aligns itself along, and follows the oscillations of, the local 
evanescent electric field vector. The appropriate Rabi frequency in this field-dipole 
orientation picture is thus given by [22] 

hnn{r) = \a\ fiEiiv) (16) 

where n and Ei are the magnitudes of these vectors. Using Eqs.(4) and (9) we can 
write the square of the Rabi frequency in the limit 7 — > in the following form, 
emphasising the dependence on the angle of incidence 0, which is related to kn through 



c^k?, = uj'^t2 sin . 



^r{^A,z) 



4|a|V'A;,22c'e-2^^i^ 



i + ^f^_iV 
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where A is a convenient scaling frequency defined by 



^1 



(1 + kJki,) (17) 



n,e^ 



A^ = ^ , (18) 

m*eoeid 

with d a convenient scahng length. Note the presence of the second term in the ex- 
pression between the last pair of brackets in Eq.(17) which is proportional to kj/k'^^. 
This represents the contribution from the z-component of Ei, and since kzi can be 
very small (when the angle of incidence is close to the total internal reflection angle), 
this could lead to an enhancement, in contrast to the contribution from the in-plane 
component (corresponding to the first term in the last pair of bracket in Eq.(17)) which 
does not exhibit this feature. 

Figure 2 displays the effects of varying the sheet electron density Ug on the value of the 
squared Rabi frequency i7|,(cj, 0, 0), i.e. evaluated at the surface z = 0. This quantity 
provides an indication of the effectiveness of the structure as an atomic mirror. The 
curves correspond to different values of angle of incidence 0. There are two features 
worthy of note here. First, for a given (p the variation of Q'j^ with increasing Ug is 
practically flat at the value corresponding to n^ = (absence of the metallic sheet) 
up to a point where it changes rapidly, increasing to a maximum and then dropping 
sharply to a relatively small value at a higher range of n^. This feature at large n^ 
signifies screening effects. The peak of ^2^ corresponds to the value of Ug satisfying 
the condition k^id = cu^/A^, which is similar to (but not the same as) the dispersion 
relation of the surface mode in this system which is [24-26] 

kzid up' 



1 ~ ^5^^2(62 -ei) 



i/2\ ^2 



(19) 



This surface mode, in turn, differs from the usual surface plasmon mode arising on a 
semi-infinite metallic half space. At an angle of incidence approximately equal to 0o for 



total internal reflection, the Rabi frequency exhibits a resonance which is at least two 
orders of magnitude larger than the value corresponding to the absence of the metallic 
sheet. At = 00 (not shown in the Fig. 2) the squared Rabi frequency Jl^ would be 
in the form of a delta function. In principle, then, pronounced repulsion effects would 
be expected in conditions corresponding to the region of the resonance. 

The real mirror action is only partially influenced by the behaviour of the Rabi fre- 
quency; dynamical effects are, of course, controlled by the forces acting on the atomic 
centre of mass. The atom experiences besides the average radiation force given by 
Eq. (p!2D , an attractive force due to interaction with the vacuum flelds which are con- 
strained by the presence of the surface. At distances from the surface large compared to 
a reduced transition wavelength Ao/27r the atom-surface force takes the Casimir- Polder 
form. At distances smaller than \q/2tt the force assumes the van der Waals form. As 
we show shortly, the important region for the atomic mirrors considered here, is, in 
fact, the van der Waals regime for which the force can be written as 

F™(^) = -^ (20) 

where t/vw is the image potential. As is well known, the van der Waals potential arises 
as the position-dependent change in the radiative self energy due to the presence of 
the surface, but it depends on the orientation of the electric dipole moment relative 
to the surface. It should be emphasised that the evanescent fleld has no other role to 
play in the determination of the van der Waals potential, except that it is responsible 
for the dipole orientation. The average electric dipole aligns itself parallel to §1, a unit 
vector in the direction of the local evanescent electric fleld so that Ei(r) = eiii^i(r). In 
the fleld-dipole orientation picture this state of affairs applies at every point along the 
atom trajectory. The leading contribution to the van der Waals potential arises from 
the interaction of the electric dipole with its image and we therefore have 

f^vw(^) = - ^/ 3 |3(ei.z)(§i.z) - ei.~ei} (21) 



where the unit vector §1 is the image of §1. 

The atomic reflection process is controlled by two separate mechanisms. Firstly, the 
spontaneous force F^ acts on the atom in the direction of ky and, secondly, the combined 
force Frf + Fvw acts to repel the atom from the surface, decelerating its motion along the 
z-axis towards the surface, attaining an instantaneous halt before changing direction 
away from the surface at an appropriate turning point in the trajectory. This behaviour 
can be seen more clearly by examining the corresponding repulsive potential Ut{"v, z) = 
Uvw{z) + f/d(v, z) which is shown in Fig. 3 for a typical set of parameters. Note that by 
controlling the average dipole direction the field-dipole orientation picture influences 
the variations of both U^„ and Ud- 

The trajectory of the atom of mass M approaching a mirror for a given set up is 
obtainable by solving the equation of motion 

M— = F, + Fd + F™ - Mgz (22) 

subject to given initial conditions. Figure 4 displays typical trajectories in the plane 
of incidence. The parameters are such that the spontaneous rate F is taken to be the 
free space value Fq. This is in fact a very good approximation in the trajectory region 
which is sufficiently far from the metallic sheet and the substrate. The static detuning 
is taken to be Aq = 5 x IO^Fq. Finally the intensity of the light is assumed to be 
/ = 2.0 X lO^W m~^. It is straightforward to verify that for this light intensity and for 
the detuning value used here, the saturation parameter S = 2r2|j(0)/(A^ + Fq) ^ 0.23 
which conforms with the low intensity regime. 

In view of Fig. 4 one concludes that the structure operates as an atomic mirror in three 
of the cases displayed, while for the fourth case the trajectory of the atom terminates 
with a collision at the surface. An approximate guide to the condition leading to a 
collision with the surface is to compare the maximum height Umax of the potential 



in Fig. 3 to the initial kinetic energy Mt>^(0)/2. For f^(0) > J2Umax/M' a colhsion 



occurs. This interpretation indeed conforms with the results of the type shown in Fig. 
3. Note that reflected atom trajectories are, in general, asymmetric with respect to the 
turning point. This is a consequence of the action of F^ which in the present example 
where the light and the atom are incident on the same side (left side) of the z-axis, 
accelerates the atom to the right along the surface. 

In conclusion, we have explored the influence of adding a metallic sheet to the usual 
evanescent mode atomic mirror set up. The theory presented provides information 
about the range of metallic sheet densities and the angle of incidence at which the 
Rabi frequency exhibits pronounced enhancement effects. We have also adopted the 
field-dipole orientation picture in which the average atomic dipole oscillates along the 
direction of the electric field of the evanescent mode at every point in the trajectory. We 
have seen that this step which allowed us to identify the direction of the average dipole 
moment vector as that of the local electric field vector at every point in the trajectory, 
has important consequences for the evaluation of the forces and, hence, the dynamics 
of the atom. Our results show that enhancement arises at an angle of incidence close 
to that of the total internal reflection condition, but requires a relatively high sheet 
density. 
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Figure Captions 
Figure 1 

Schematic arrangement of the elements comprising an evanescent mode atomic mirror 
with a metalhc sheet. The plane of incidence contains the internally reflected light 
beam as well as a typical atomic trajectory. Here both the atom and the light are 
initially propagating on the same side of the vertical (z-axis). 

Figure 2 

Variation with areal density Ug (in units of nf^'"^'^ = 5.573 x lO^^m"^) of Q'j^, the squared 
evanescent mode Rabi frequency (in arbitrary units), evaluated just outside the metallic 
sheet in the limit of small 7. The different curves correspond to different angles of 
incidence = 42.00° (dashes); = 47.00° (dash-dots); = 52.00° (dash-double dots). 
The inset (solid curve) shows the variation for = 41.37°, which is very close to 
the angle 0o for total internal reflection. The parameters are such that 62 = 2.298, 
corresponding to an angle of total internal reflection 0o = 41.27°; Fq = 6.128 x 10® s~^; 
Ao = 5.0 X lO^Fo and the intensity of the light is taken as / = 2.0 x lO'^Wm"^ for the 
Rb resonance at the transition wavelength Aq = 780nm. The scaling length d is taken 
as (i = 15 nm. 

Figure 3 

Variation with z (in units oi d = 15nm) of the combined static dipole potential U = 
Ud + Uvw acting on a Rb atom (in units of hTQ/2). The curves correspond to a fixed 
value of ris = 215nf '"'^'', but different values of angle of incidence 0. They are as follows: 
= 42.00° (dashes); = 45.00° (dash-dots); = 47.00° (dash-double dots); = 41.37° 
(solid curve). The dotted curve shows the variation of the van der Waals potential. In 
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the evaluation of the potentials, the direction of the dipole moment vector conforms 
with the field-dipole orientation picture. The parameters are the same as those used 
in the evaluation of Fig. 2. 

Figure 4 

Atomic trajectories of a Rb atom in the atomic mirror arrangement shown in Fig.l with 
the x-z plane as the plane of incidence. The angle of incidence is fixed at = 42.00° and 
the metallic sheet corresponds to n^ = 215nf^'"'^^, as in Fig. 3. In all cases the initial 
position of the atom is at the point {x = 0,z = lOOd), where d = 15nm. The different 
trajectories correspond to the same initial condition for the horizontal component of 
the velocity f||(0) = 0.054 ms^^, but differ in their initial z-component of velocity ^2(0). 
They are as follows: fz(0) = O.lms"-^ (solid curve); ^2(0) = 0.4ms~"'^ (dashed curve); 
■^2(0) = O.Tms^"*^ (dash-dotted curve); ^2(0) = l.lms^^ (dash-double dotted curve). 
The parameters are the same as those in Figs. 2 and 3. 
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